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We report on a high statistics lattice QCD calculation of the mass of the rj and rj' mesons using 
ASQTAD improved staggered fermions. The calculation used two ensembles with different lattice 
spacings and pion masses. We also report results for rj-rj' mixing. The results are in satisfactory 
agreement with other lattice calculations using other fermion formulations and with experiment, 
given the the unphysical quark masses used. We see no evidence of abnormal behavior at the lattice 
spacings studied. 



I. INTRODUCTION 

The explanation for the large mass of the rj' meson 
relative to the masses of the light pseudo-scalar mesons 
involves nonperturbative physics, such as topology and 
the J7(l) anomaly so these masses present an in- 

teresting challenge for lattice QCD. There are many phe- 
nomenological applications that require knowledge of the 
properties of the 77 and 77' mesons, such as their decay con- 
stants. For example, Di Donato et al. [3] have recently 
noted that a theoretical description of 77-77' mixing is re- 
quired for understanding the CP asymmetries of charmed 
and bottom mesons with 77 or 77' in the final state. Fleis- 
cher et al [si have recently studied the dependence of 77-77' 
mixing on the decay ^ — > J /ipr]'^'^ , and which could 
be used to look for physics beyond the standard model. 
There is currently an active experimental program ^] 
into the properties of the 77 or 77' mesons at experiments 
such as WAS A ^, KLOE-2 and MAMI [9^. Lat- 
tice QCD can in principle compute the 77-77' mixing from 
first principles, so it should be able to help interpret the 
experiments. 

The improved staggered fermion formulation, known 
as ASQTAD fermions [13, [HI, has resulted in some very 
accurate lattice QCD calculations which have been val- 
idated a gain st experiment with a high degree of preci- 
sion [1^ Il3j|. The results from the improved staggered 
fermions physics program have been recently been re- 
viewed by the MILC collaboration 14 1. 

Despite the success of the phenomenology from stag- 
gered fermions, concerns have been expressed over the 
theoretical basis of the formalism. The central issue 
arises from the so-called 'rooting' of the staggered de- 
terminant This is invoked in order to correct 
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the number of fiavors/tastes of sea quarks which would 
otherwise arise. These aspects of the staggered fermions 
formalism have been reviewed by Diirr [17] , Sharpe [l8| , 
Creutz 19] and Kronfeld 

The physics of the rj' meson is sensitive to the topol- 
ogy of QCD. Since some of the concerns about rootin g th e 
staggered formalism are related to the topology [TGlflQl. 
the rj' could be a place where incorrect results are ob- 
tained. The properties of the rj' meson are also thought 
to be related to the eigenvalue spectrum of the Dirac op- 
erator. The eigenvalue spectrum of a type of improved 
staggered fermions called HISQ has recentl y b een shown 
to agree with continuum expectations [2l| - |24| . which is 
suggestive that the physics of the 77' meson is correct with 
improved staggered fermions. The recent paper by Don- 
ald et al [23] also contains a rebuttal of some previously 
raised theoretical concerns. A calculation of the proper- 
ties of flavor singlet mesons is a further and crucial test 
of the validity of the improved staggered formalism. 

The very successful physics program that uses im- 
proved staggered fermions has largely been independent 
of the properties of the 77 and 77' mesons [3] • The MILC 
collaboration has recently used the properties of the 77' in 
their study of the topological susceptibility [1^. The 77 
and 77' mesons could indirectly affect other lattice calcu- 
lations. For example, the 77 and 77' mesons are common 
decay products of hadrons that decay via the strong force, 
such as the an meson [26l . [27| or the proposed exotic hy- 
brid mesons [23]- However the majority of lattice QCD 
calculations performed by the MILC collaboration are 
probably independent of the properties of the rj and 77' 
mesons. 

In [11] we reported on methods to compute discon- 
nected diagrams with improved staggered fermions. In 
that study we found that we needed much higher statis- 
tics than those of a typical ensemble generated by the 
MILC collaboration (at that time). This new work is 
based on generating many more configurations in order 
to improve the statistics. We have recently reported [13] 
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results for the masses of the 2++, 0++ and gluebahs 
obtained as a byproduct of this calculation. 

Lattice QCD calculations of flavor singlet quantities 
are not as well developed as those of flavor nonsin- 
glet quantities, because they are computationally more 
expensive. There have been lattice QCD calculations 
of the 77' mesons with 2+1 flavors of sea quarks, by 
the JLQCD collaboration [3l|, RBC/UKQCD collabo- 
ration [ill, Hadron Spectrum CoUaljoration [s^, and 
TWQCD and JLQCD cohaborations [p]. A preliminary 
study of the form factors for the semileptonic decay of 
Dg -> rjg has been carried out [35|. There is also a pre- 
liminary calculation of the masses of 77 and rj' mesons, 
with 2+1+1 flavors of sea quarks, from the ETM collab- 
oration [36| . 

There have been a few lattice QCD calculations of the 
flavor singlet pseudoscalar meson with n/=2 sea quarks. 
However, the strange quarks in the sea play an essential 
role in lattice QCD calculations of the mass of the rj' 
meson. For example Jansen et al. (STj found the ground 
state of the flavor singlet pseudoscalar meson in = 2 to 
be 0.865 (65) (65) GeV compared with the experimental 
value of the r] mass of 0.548 GeV which lattice QCD 
calculations with 2+1 sea quarks should reproduce. 
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FIG. 1. Normalized histograms of raw disconnected corre- 
lator measurements with local (left) and fuzzed (right) pseu- 
doscalar source operators. Measurements are for the light- 
light correlator on the coarse ensemble. Note that symme- 
try increases with At and the most likely value remains at 
D{At) = for all At. Histograms contain 335168 measure- 
ments (5237 configs x 64 timeslices. The heavy dashed green 
curve is a fit of Eq. ^ to the fuzzed At = 10 histogram.). 



II. DETAILS OF THE LATTICE 
CALCULATIONS 



We used the improved staggered fermion action called 
ASQTAD [13, im with the tadpole improved Symanzik 
gauge action. We generated two ensembles at two lat- 
tice spacings, using QCDOC machines [38|. The basic 
parameters of the lattice QCD calculations are in Ta- 
ble ID We generated wfj"*''''' configurations, but base our 



final calculation on a subset of N^^'^'^^^ configurations. 
The RHMC algorithm was used to generate the configu- 
rations [H Ei]. The details of the tuning of the RHMC 
algorithm are in [4l| . 

The use of new measurement techniques has improved 
the precision of many lattice QCD measurements. How- 
ever, for many applications of lattice QCD the use of 
higher statistics is still an essential requirement. In par- 
ticular, this is so for glueball studies [4l], or for anything 
which includes the calculation of disconnected diagrams. 

Our aim was to generate 30,000 trajectories. For com- 
parison the ETM collaboration has typically used ensem- 
ble sizes of 5000 or 10000 trajectories [s^- The length of 
simulation required also depends on the autocorrelation 
time. There have been some recent high statistics lattice 
QCD calculations using an anisotropic lattice (43j . 

In our methods paper [2^ on the 77' we found that 
the correlators did not have a Gaussian distribution. A 
similar observation has been made by the ETM collabo- 
ration in their calculation with n/ = 2 flavors of sea 
quarks. In this calculation we again observe a long tail 
in the distribution of the disconnected correlators. 



In [29| we pointed out that the long-tailed distribu- 
tion is an inherent challenge to measuring disconnected 
correlators. Under the assumption (well-justified by nu- 
merical observation) that the pseudoscalar fermion loop 
operators fluctuate within a Gaussian distribution about 
zero, the product of two loop operators, measurements 
disconnected correlator d must fall in a long-tailed dis- 
tribution, specified by 



n{d) = Aejipi~Bd)Ko{C\ 



(1) 



where Kq is a modified Bessel function of the second kind, 
and A, B and C are constants. 

The distribution of disconnected correlator measure- 
ments will always have its peak at zero, with the signal 
coming from the asymmetry of the long tails. This be- 
havior may, however, be masked by any binning of mea- 
surements. We show some histograms of disconnected 
correlator measurements in Fig. [TJ 

The UKQCD/RBC collaboration has, however, not ob- 
served any "large statistical excursions" in the distribu- 
tion of the disconnected correlators in an unquenched 
calculation with 2+1 fiavors of sea quarks [32| . 

It might therefore be tempting to suggest that the long 
tail in the distribution of the disconnected correlators is 
caused by the lack of 'lattice' chiral symmetry in the 
staggered and twisted mass fermion operators, compared 
to the domain wall formalism. However, we explicitly 
checked that we used sufficient numbers of noise sources 
so that our errors were dominated by the gauge noise. 
The wall sources used by the RBC/UKQCD collabora- 
tion [32] can be thought of as using a single noise source 
but with partitioning in time. It might be that the dis- 
tribution of the errors on disconnected correlators in the 
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Name 




uo 


xT 


ami 


arus 


^(total) 
cfK 


iVtraj 


^(mcas) 


ro/a 


coarse 
fine 


6.75 
7.095 


0.8675 
0.8784 


24^ X 64 
32=* X 64 


0.006 
0.00775 


0.03 
0.031 


5237 
2867 


31422 
17202 


4452 
2808 


3.8122(74) 
5.059(10) 



TABLE I. Summary of ensembles used in this calculation. We use the same convention for the quark masses as used by the 
MILC collaboration. 



latter case are not dominated by the gauge noise. This 
issue is discussed by Wilco x and there is a recent 
study by Alexandrou et al. [45| . 

Recently Endres et al. [46] have discussed heavy-tailed 
correlator distributions in condensed matter simulations 
and lattice QCD calculations ^^'^ started to develop 
new techniques to extract physics in such situations. 

We note that the width of the disconnected correlator 
distribution increases with fuzzing of the sources. The 
taste-singlet pseudoscalar operator in the staggered for- 
mulation has the quark and anti-quark separated by four 
links. It is possible that the introduction of the gauge 
links in the covariant shift or in the fuzzing increases the 
width of the Gaussian distribution of the loop operators 
and hence the width of distribution of the disconnected 
correlator described by ([T]). In other formulations, the 
distributions may possibly be narrower and less trouble- 
some. 

It can be difficult to make a good choice of strange 
quark mass to use in a lattice QCD calculation, although 
the choice is getting easier with increased experience of 
2-1-1 calculations. In the original calculations by the 
MILC collaboration, the estimate of the strange quark 
mass turned out to be off by as much as 25% [lj|. The 
MILC collaboration then corrected for the mismatch of 
the strange quark mass by using their extensive partially 
quenched data sets and later unquenched runs. Their 
updated estimates for the strange quark mass (in lattice 
units) were 0.035(7) and 0.0261(7) on the coarse and fine 
ensembles respectively 14]. However we have used the 
new value of mg that MILC used on the coarse [1^, so 
as to interpolate with their original value of the strange 
quark mass. Indeed, our coarse ensemble has already 
been used in an analysis of the strangeness content of 
the nucleon, by the MILC collaboration (49| . 

We first discuss the interpolating operators and two- 
point functions for flavor singlet pseudoscalar mesons us- 
ing Wilson fermions with n/ degenerate quarks. A longer 
discussion can be found in our original paper [29| . Two- 
point functions 



G{x',x) = (0(.T')0(a;)) 
are constructed using the interpolating operators 



0{x) = ^i/',(a;)75?/'j(a;) . 



(2) 



(3) 
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FIG. 2. Effective mass plot for the connected part of (75 x 1) 
operators for the fine ensemble. 



(D) part 



G{x' , x) = nfC{x' , x) — n'jD{x' , x) 



(4) 



Using standard Wick contractions, the two-point func- 
tion decomposes into a connected (C) and disconnected 



In the staggered formalism the native four tastes of sea 
quarks are reduced to n degenerate flavors through the 
j-rooting of the fermion matrix. The four valence tastes 
manifest themselves in valence loops, which come in too 
large by a factor of four. A connected correlator is es- 
sentially a single valence loop and requires normalizing 
by a factor of j to get the single-flavor contribution. A 
disconnected correlator contains two valence loops (along 
with an arbitrary number of sea quark loops). Comput- 
ing the single-flavor disconnected contribution therefore 
requires normalizing by two factors of j. This factor is 
sometimes referred to as 'valence rooting' [l^. We in- 
clude these factors implicitly in our definition of D and 
C. 

In the stag gered formulation we use the Kluberg-Stern 
notation . The Goldstone pion operators are (75 x 75 ) , 
and the flavor singlet pion operators are (75 x 1), where 
the first index refers to spin and the second index to taste. 

The signal-to-noise ratio for the flavor singlet pseu- 
doscalar mesons rapidly falls with time, so it is impor- 
tant to use spatially smeared operators to project onto 
the ground and first excited state as soon as possible. 
We used the fuzzing technique [sij with a fuzzing length 
of 4 (this must be even to respect the symmetries of the 
staggered action). 
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lattice spacing we find the value of the Goldstone pion 
mass on the coarse and fine ensemble to be 295 and 349 
MeV respectively. 

The mass of the connected strange-strange pseu- 
doscalar meson rrirj^ = 0.6858(40) GeV, has recently been 
determined by the HPQCD collaboration 56]. By using 
m^^ cx TOs we find that the mass of the strange quark is 
mistuned by -20% on the coarse ensemble and -1-5% on 
the fine ensemble. 

The masses in table |TT] can be used to check the re- 
duction of taste breaking as the lattice spacing is re- 
duced [2^. Defining 



d= (m^ 



(75x1) 



''(75x75) 



(9) 



FIG. 3. Effective mass plot for the connected part of (75 x 1) 
operators for the coarse ensemble. 



we find dfin^Aicoarse = 0.39(1). The MILC collaboration 
obtained [20| dfine/dcoarsc — 0.38(3) at similar param- 
eters. Note that MILC use ri in equation [9l but the 
conversion between ri and tq drops out in the ratio. 



We computed a variational correlator matrix C using 
four basis states - light and strange quark pseudoscalar 
correlators combined with local and fuzzed sinks and 
sources: 

C(% = (0.(00, (0)t). (5) 

In Figs. [2] and |3] we show the effective mass plots for 
the local and fuzzed smeared operators for the (75 x 1) 
pion operator. For the main rj-rj' analysis, we use the 
factorizing fit model (52j shown in Eq. 

N 
k—l 

where i,j run from 1 to 4. Here N is finite and all states 
are treated as stable. In final fits we used no Bayesian 
constraints. 

Physics can also be extracted from the measured 
variational matrix using the generalized eigenvalue 
method [il-Ill. 

C{t)v„{t, to) = A„(i, to)Cito)Vn{t, to) (7) 

where w„(i,to) sue the generalized eigenvectors. Typi- 
cally a small time to is chosen. 

The energies are extracted from the generalized eigen- 
values 

A„(t,to) = e-^"(*-*°). (8) 

We use the effective mass from Eq. ([5]) as a check on the 
masses from the fits to Eq. ([5]). 

We used the recently updated value of ro=0. 4661(38) 
fm [5^ to determine the lattice spacing from our mea- 
surements of ro/a [3^ . The inverse lattice spacing deter- 
mined from ro is 1.61 and 2.14 GeV for the coarse and 
fine ensembles respectively. With our estimates of the 



III. ANALYSIS OF DISCONNECTED LOOPS 

In [29I we reported on methods to compute dis- 
connected diagrams with improved staggered fermions. 
We found that the technique proposed by Kilcup and 
Venkataraman [53] was the most efficient one of those 
we tested for computing the disconnected diagrams of 
(75 X 1) and (1x1) operators. We therefore used that 
with 32 Gaussian noise vectors in this calculation. 

The integrated autocorrelation time for the light and 
strange (75 x 1) loops are 42(7) and 24(3) in units of 
trajectories respectively. This can be compared with the 
autocorrelation times of the connected correlators for Eu- 
clidean time separation 3 of 14(4) and 12(5) for the light 
and strange quarks respectively. 

A. Ratio of disconnected to connected diagrams. 

An important test of the staggered formalism is to 
check the long time behavior of the D{t) /C{t) ratio of the 
correlators in Eq. (g]). In quenched QCD, the D{t)/C{t) 
ratio rises linearly with time but, in unquenched QCD, 
the ratio should tend to 1 for large times. This test was 
performed for clover fermions in (58| . In our earlier study 
(29| ] , the statistics were not large enough to determine the 
large time behavior of the ratio in unquenched QCD. In 
Fig. m we plot the D{t)/C{t) ratio for the coarse and fine 
ensembles respectively, as well as a quenched ensemble 
with the same lattice spacing as the coarse ensemble. 

With the higher statistics afforded by the two new 
ensembles there is now a clearer difference between 
the quenched and unquenched D{t)/C{t) ratios. The 
D{t) /C{t) ratio for the two unquenched seem to be reach- 
ing a plateau close to 1. As explained in Sec. HIl the 
staggered ratio depends on the 'valence rooting' factors 
of 1/4, so this is a good test of the staggered formalism. 
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nge 
(75 X 1) 


6.75 
7.095 


0.183(1) 
0.1632(7) 


0.325(2) 
0.206(1) 


0.387(2) 
0.3283(7) 


0.4683(8) 
0.3505(6) 



TABLE II. Masses for pseudoscalar mesons in lattice units from the connected correlators. 
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FIG. 4. D(t)/C{t) ratio for the coarse and fine ensembles. 
We also include the result from a quenched lattice QCD calcu- 
lation [29j ] with a lattice spacing close to the coarse ensemble. 



FIG. 6. Comparison of fit results (horizontal bands) to 
effective masses (with At — t — to — 2) from the variational 
method for a 4 x 4 correlator matrix on the fine ensemble. 
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FIG. 5. Comparison of fit results (horizontal bands) to 
effective masses (with At = t — to — 2) from the variational 
method for a 4 x 4 correlator matrix on the coarse ensemble. 



B. Fit results for the masses 

We plot the fit results against the effective masses from 
the variational method ([5]) in Figs. [S] and IHl Here, we 
set A = i — io = 2in(l8]). The parameters from the final 
factorising fits are in given Table Hill In arriving at these, 
we varied tmin and iVoxp looking for a sweet spot of sta- 
bility with respect to these, for x^/dof < 1, and for error 
bars small enough to distinguish three states. Plateaux 
in tmin were small due to the competition between ex- 
cited states at small t and the onset of large noise/signal 
at larger t. An SVD cut of around 10 ~^ was necessary 



to get fits with reasonable confidence levels. 

Note that the statistical errors on the final mass pa- 
rameters deduced from these correlated fits are generally 
somewhat less than the individual variational effective 
masses calculated at each value of t, even though the 
same number of configurations were used. The global 
factorizing fits include all the data from the different 
flavor and fuzzing channels simultaneously, and from a 
large range of i-values. They also take correlations into 
account. The individual effective mass estimates from 
([5]), on the other hand, are each determined by relatively 
small samples of this data set — uncorrelated measure- 
ments from two time-slices yield independent estimates 
for each eigenvalue. 

We found it necessary to restrict imax to moderate val- 
ues (16 and 20) since, at larger values, the correlators 
can fluctuate below zero. We do not attribute this large-t 
behavior to unphysical effects. Rather, the unique sta- 
tistical properties of disconnected correlators mean that 
a small number of configurations at the tail of the dis- 
tribution can cause large fluctuations of the mean, as 
discussed in [29| . 

We associate the lowest state with -q and the first ex- 
cited state with the ry'. 

IV. COMPARISON WITH EXPERIMENT AND 
OTHER LATTICE CALCULATIONS 

Given that we only have two masses at two different 
lattice spacings, we cannot make a controlled continuum 
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TABLE III. Masses for taste singlet pseudoscalar mesons in lattice units. The errors are from the jackknife method. The 
corresponding fit methods are described in the text. 
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FIG. 7. Summary of mass dependence of the rj and rj' as 
a function of the square of the pion mass. We also include 
data from the UKQCD/RBC [32l, Hadron Spectrum Collab- 
oration [s^]. The state listed as 77(3rd) is the next excited 
state in our fits. 
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FIG. 8. Summary of the rj and ry' as a function of the 
square of the lattice spacing. We also include the mass of the 
connected (75 x 1) operators to show the scale of the lattice 
spacing errors. 



or chiral extrapolation. However, we do wish to check 
that the resulting parameters are physically reasonable. 

In lattice QCD calculations with two flavors of sea 
quarks, it was found that the mass of the 772 light flavor 
singlet pseudoscalar meson had very little mass depen- 
dence |37| . 

In Fig. [7] we plot our results for the masses of the 
77 and rj' mesons as a function of the square of the 
pion mass. We also include the recent results from the 



UKQCD/RBC collaboration [321 ^^'^ the hadron spec- 
trum collaboration [33| . Note that while the lattice spac- 
ing is larger for the coarse ensemble than the fine en- 
sembles, the light quark mass is smaller for the coarse 
ensemble than the fine. 

The mass from the (75,1) pion operator is different 
from the Goldstone (75,75) pion operator, because of 
taste breaking. In Fig. [8] we plot, as a function of the 
square of the lattice spacing, the masses of the and rj' 
mesons along with the mass splitting between the (75, 1) 
and (75,75) pions. 
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MIXING 



In a quantum mechanics analysis of the mixing of the 
7] and rj' [sol . [60j . the starting point is two bare flavor 
singlet pseudoscalar mesons. As interactions are intro- 
duced, the two flavor singlet mesons mix and a mixing 
angle is introduced. This mixing angle is a parameter in 
the various decays of the rj and rj' mesons. 

In unquenched lattice QCD calculations it is not clear 
how to start from unmixed states. It is possible to esti- 
mate the elements of the mixing matrix in quenched and 
partially quenched QCD [6l|, H^j. We also note that the 
properties of the connected strange-strange pseudoscalar 
meson determined by the HPQCD collaboration [56'| are 
also useful for developing a partially quenched mixing 
analysis [6]| . 

In QCD, the mixing between the 77 and rj' mesons is 
understood at the amplitude level. In phenomenological 
analysis of decays, the mixing is parametrized in terms of 
decay constants and 1 or 2 angles. Many of the analyses 
of the decays also include additional strong interaction 
physics, such as form factors [63| . 

Di Donato et al. [1] and [sj Fleischer et al. have re- 
cently discussed the importance of the mixing of 77 and 
tY mesons on various decays. Shore and collaborators 
[3] have derived equations relating the masses and decay 
constants of the 7] and ry' to the topological susceptibility. 

The phenomenological analysis of 77-77' mixing is usu- 
ally described in terms of the axial decay constants [l|. 
To study the decay constants of flavor singlet meson oper- 
ators with staggered fermions requires the (7475, 1) pion 
operators. These were not included in our calculation. 
Since there is no Ward identity for the (7475, 1) operators, 
renormalization factors would need to be computed [i^l . 
The connection between i] — rj' mixing with axial cur- 
rents or pseudoscalar currents is discussed in section 3.8 
of Feldmann's review [1]. To first order the mixing angle 
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is the same with axial or pseudoscalar currents. 

There are two popular bases in which to express rj/rj' 
states. One basis is the S'/7(3) basis: 



1 



1%) = ^ + \dd) + \ss)) 

hs) = ^ {\uu) + \dd) - 2\s-s)) 



(10) 



The rj and rj' mesons states are then linear combinations 
of the basis states: 



cos Op —smOp \ I I ryg) 
sin 9p cos Op 1 \ \ 770) 



(11) 



The r] and rj' mesons can also be described using the 
quark flavor basis: 



1 



\Vq) - + \dd)) 



\Vs) = \ss) 




cos (pp — sm (pp 
sin (f) p cos 0p 




(12) 



(13) 



The two mixing angles are related to each other: 

Op^ffyp- arctan ^/2 ^ (pp - 54.7° . (14) 

In lattice QCD calculations it is not possible to start 
from pure basis states such as \'qq) and |?7s). Instead the 
mixing is described in terms of amplitudes. For example 

via 



we can define an amplitude a^^^ 

a^jg^j = (0 I ?78 I ri). 



(15) 



Other amplitudes, for different interpolating operators 
and mesons, are similarly defined. For example, the am- 
plitudes ^ and a,,^^ are parameters in the fit model in 
Eq. ([6]) , with numerical indices replaced by the names of 
the meson or interpolating operator. 

With two possible interpolating operators and two 
states, there are in principle four possible decay con- 
stants. The discussion of mixing in quantum mechanics 
suggests that these four amplitudes can be parameterized 
in terms of two decay constants and a mixing angle. 

It has been found that, with S'?7(3)-breaking, there 
cannot be one single mixing angle in the SU (3) basis |6al. 
For example from a review of the literature, Feldmann [1| 
notes a spread of 10 degrees in the mixing angles ex- 
tracted from experiment. In the SU{Z) flavor basis the 
mixing of decay constants is expressed by 



air) ao,, \ / /s cos 6*8 -/o sin 6*0 
as,,' aor,' / I /s sin 6*8 /o cos 6*0 , 



(16) 



The decay constants can also be computed in the quark 
basis: 





ClsTj \ 




\aqri' 







fq cos (f)q -/sSinc 

fq sin (pq fs COS (j) 



(17) 



To extract the mixing angles from the matrix in 
Eq. p7)) . the following combinations can be used: 



tan 



tan ( 



fq' 



(18) 



(19) 



There are arguments [U,!!^ that suggest (f>q ~ (j)s- Feld- 
mann [1| reviews the determination of the mixing angles 
from various processes and he finds | <l)q ~ cf)s |< 5 de- 
grees. The UKQCD/RBC colaborations use perturba- 
tion theory to identify conditions where there is only one 
angle [32j . 

In Fig. [5] we show the fit amplitudes plotted in couplets 
{cLri , ttrf' ) demonstrating the mixing angles in the quark 
flavor basis as in Eq. (IT71) . 

Because of SU(3) symmetry breaking, fq is not equal 
to fs in Eq. ([TT]), hence a rotation of the interpolating 
basis from the quark to the flavor basis will no longer give 
a matrix of the same form as Eq. (|17p . The use of mixing 
angles to parametrize the mixing matrix is similar to the 
polar decomposition of a matrix j67| , where a matrix can 
be written as the product of a Hermitian matrix and a 
unitary matrix. In this application the Hermitian matrix 
is actually diagonal. 

The RBC/UKQCD and Hadron Spectrum Collabora- 
tion extract a single angle by computing the following 
angle 



tan 



2 I est 



(^qrjO'srj' 



(20) 



so that tancff^^ is the geometric mean of tant/)^*** and 
tan ipf^ . 

The relation (pq sa (j)^ is an assumption, under which 
the four possible parameters are reduced to three, hence 
it is interesting to test it using lattice QCD. We have de- 
termined the mixing angles using different methods and 
present the results in Table |lVl The angles (jf^^, (j)°^^ and 
(ff^^ were obtained from the fits to Eq. ([6]), using Eqs.([T8|. 

and (PO)) . Local operators were used for this. We 
also found similar results using fuzzed operators. 

The angle (j)^\ in Table |TVl was determined by simul- 
taneously fitting the smeared and local matrix of ampli- 
tudes in Eq. P7|) using a single angle 0^' = 0^ = 0s, 
with two decay constants for the local matrix, and two 
decay constants for the fuzzed correlators. 

One cannot confidently conclude from the results in 
Table [TVl that a single mixing angle does not describe the 
data. As an exercise we ignored the statistical errors on 
the amplitudes from the /3 = 7.095 ensemble and rotated 
the matrix in Eq. [17] by an angle a in the interpolating 
operator basis and computed (ff^^*" and (fif*'. Around 9 
degrees from the quark basis we found (jf^^ = 0°*** as a 
function of a. This shows as expected that the quark 
basis is a good basis for having a single mixing angle. 
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/est 


/est 


/est 

<? 






6.75 


25(4) 


36(2) 


31(4) 


34(3) 


8.2/3 


7.095 


40(5) 


34(2) 


37(3) 


34(2) 


3.7/3 



TABLE IV. 77-77' mixing angles in degrees defined and determined as described in the text. 







coarse local 


q 




A 


coarse local 






© 


coarse fuzzed 


q 






coarse fuzzed 









fine local 


q 




□ 


fine local 









fine fuzzed 


q 




a 


fine fuzzed 


s 






(|) ^ 34 deg 







-0.2 -0.15 







quadratic GMO 
radiative decays 

radiative decays + glueball mixing 

charm-eta production 

photon fusion 

UKQCD/RBC 

HSC 

UKQCD (2000) 
UKQCD2011 (this work) 



0.4 



0.5 



(GeV) 



FIG. 9. The fit amplitudes plotted in couplets, (a^,a^/), with 
the rj amplitude as the horizontal coordinate and the 77' ampli- 
tude as the vertical component, is a graphical representation 
of Eq. p7p . This illustrates several nontrivial results, includ- 
ing the consistency of the mixing angles from local (open) 
and fuzzed (filled) sources, and the consistency of mixing an- 
gles from coarse and fine simulations. The extent to which 
the light quark and strange quark branches lie at right angles 
to each other is a measure of the aptness of the single mixing 
angle description. We use lattice units; the inclusion of renor- 
malization factors will not change the angular distribution. 



This method can be used with future more precise lat- 
tice QCD calculations, to determine the optimal basis of 
interpolating operators, such that the amplitudes can be 
described by a single mixing angle. 

In Fig. [TU]we plot the mixing angle in the quark basis 
((^^'), obtained using the fit method, as a function of the 
square of the pion mass. We also include some experi- 
mental numbers for the mixing angle from the summary 
in [3| . We see that our results from staggered fermions are 
qualitatively in agreement with the results of the other 
two lattice groups and experiment. 



VI. CONCLUSIONS 

We have reported a lattice QCD calculation of the 
masses and mixing of the 77-77' mesons. Our key results for 
the masses (Fig.|5]) and mixing angles (Fig. [TU|) . and the 
ratio plot (Fig. 2] ) do not give any grounds for concern 
over the validity of the staggered fermion formulation. 
Our results are qualitatively in agreement with experi- 
ment and with lattice analyses using other formulations. 

Using a technique specific to staggered fermions [s^l 
we were able to compute the disconnected diagrams at a 



FIG. 10. A summary of 77-77' mixing angles as a function 
of the square of pion mass. We also include data from the 
UKQCD/RBC [s^l, Hadron Spectrum Collaboration [s^l, and 
UKQCD collaborations The UKQCD/RBC number is 

quoted in the SU(3) basis, so we have used equation [14] to 
convert it to the quark basis. 

cost roughly 10 times the cost of the connected correla- 
tors. However, as always, higher statistics are required 
for lattice QCD calculations which require disconnected 
diagrams. 

The main shortcoming of the calculation is that we 
have used only one quark mass at each of two different 
lattice spacings. It is clearly important to extrapolate 
the masses of the 77 and 77' mesons to the physical quark 
masses and to take the continuum limit. To do this, it 
may be better to use a staggered fermion formalism with 
reduced flavor symmetry breaking. The MILC collabo- 
ration 68] and HOT collaboration [6^ have started to 
simulate the HISQ 70] fermion action, which uses an 
additional level of fat links over the ASQTAD action so 
reducing the mass splitting between the masses of the 
(75 X 1) and (75 X 75) pions. 

The stout link staggered action used by Aoki et 
al. [71I [ll] could also be used, since this formalism has 
reduced taste splitting compared to ASQTAD and the 
quark masses are within 10% of the physical light and 
strange masses. 

For the phenomenology of 77 and 77' it is important to 
study the decay constants and mixing angle of the 77 and 
77' mesons. Although staggered fermions are computa- 
tionally cheap and have a Ward identity for the Gold- 
stone pion operator, it is not clear whether we can obtain 
a strong enough signal for the (7475 x 1) operator to be 
able to extract decay constants [64] . 

Even more challenging would be to study the mix- 
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ing of the the Vj and t] with the pseudoscalar gluebalL 
This requires an extension of the formahsm for mixing 
described in section |Vl In principle, the mixing angle be- 
tween the 77 and 77' and the pseudoscalar glueball can be 
obtained by fitting various branching ratios. However, 
a clear picture from the phenomenological approach has 
yet to emerge d, H, [tI, [tJ, so a first principle calcu- 
lation from lattice QCD would be valuable. The KLOE 
experiment [sj has analysed their branching fraction data 
for vector meson radiative decay to pseudoscalar mesons 
with a model and obtained an estimate of jy'-glueball mix- 
ing which is three sigma from zero. The KLOE-2 ex- 
periment plans to reduce the experimental errors on the 
analysis [1]. Given that we have already published [s^ 
a study of the pseudo-scalar glueball on these ensem- 
bles, we could at the very least include the gluonic and 
fermionic pseudoscalar interpolating operators in a single 
variational calculation (as was done for scalar operators 
in [Zl). 

This calculation would be interesting from a nonper- 
turbative QCD perspective, but also would be important 
for exploring CP violation in the decays of B and Bg 
mesons to final states that include the r\ and 77'. 

Lattice QCD calculations of the t] and 77' mesons in- 
volve an interesting mixture of: toplogy, anamolies, mix- 



ing via quark loops, and have important applications to 
phenomenology. The relatively high cost of lattice QCD 
calculations of the 77 and 77' mesons has meant that they 
are not as highly developed as the calculation of flavor 
nonsinglet hadrons pUt]. Since the computational costs 
are becoming less exorbitant, we expect that lattice cal- 
culations of the 77 and 77' mesons may soon become pre- 
cision studies. 
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